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We study the dynamics of current-biased Josephson-junction arrays with a magnetic penetration 
depth Ax smaller than the lattice spacing. We compare the dynamics imaged by low-temperature 
scanning electron microscopy to the vortex dynamics obtained from model calculations based on 
the resistively-shunted junction model, in combination with Maxwell's equations. We find three 
bias current regions with fundamentally different array dynamics. The first region is the subcritical 
region, i.e. below the array critical current I c . The second, for currents / above I c , is a "vortex 
region", in which the response is determined by the vortex degrees of freedom. In this region, the 
dynamics is characterized by spatial domains where vortices and antivortices move across the array 
in opposite directions in adjacent rows and by transverse voltage fluctuations. In the third, for 
still higher currents, the dynamics is dominated by coherent-phase motion, and the current-voltage 
characteristics are linear. 
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I. INTRODUCTION 

Studies of Josephson-junction arrays are of interest to 
model vortex dynamics as well as for their application in 
superconducting electronics Most experimental re- 
sults dealing with (vortex) dynamics in 2D arrays were 
obtained by measuring the current-voltage (I-V) charac- 
teristics. Such measurements do not give the spatially re- 
solved information needed for an unambiguous determi- 
nation of the detailed microscopic dynamics underlying 
the measured response. By contrast the low-temperature 
scanning electron microscopy (LTSEM) is a technique 
that allows for spatially resolved investigation of the dy- 
namical states in superconducting systems. As such it 
offers the possibility to determine the microscopic nature 
of the dynamics, if one correlates the experimental in- 
formation with the microscopic dynamics obtained from 
model calculations. 

Here we present both experimental images and the- 
oretical vortex dynamics results in dc-biased 2D classi- 
cal arrays (i.e. with a Josephson coupling energy Ej 
much larger than the charging energy E c ) with high 
damping and a small magnetic penetration depth, and 
in zero applied magnetic field. We gain detailed in- 
sight in the spatially resolved vortex dynamics by supple- 
menting and comparing the LTSEM results to those ob- 
tained from model calculations employing the resistively- 



shunted junction equations together with Maxwell's 
equations. Due to the small magnetic penetration depth 
of our samples, the applied dc current gives rise to 
strong induced magnetic fields at the edges of the arrays. 
These fields in turn facilitate the penetration of (current- 
induced) vortices at the edges. In this respect there is 
some correspondence between inductive overdamped ar- 
rays and continuous superconducting thin films. Both 
systems show a current-induced resistive state due to the 
nucleation of vortices of opposite vorticity at opposite 
ends of the sample and subsequent vortex motion into 
the sample, as is described for superconducting bridges 
in References [^]||- 

We include the mutual inductances between array cells 
in order to take into account the self-induced fields in 
the model simulation. We use an algorithm developed 
recently in Ref. 0], that takes into account an approx- 
imate full-range inductance matrix. A particularly in- 
teresting region of array dynamics is the current region 
slightly above the array critical current. There we find 
an intricate structure in the I-V curves (5j-|). Recent re- 
sults of LTSEM experiments in this current region were 
interpreted in terms of the collective motion of current- 
induced vortices JBj. Subsequently the detailed form of 
this motion was found to be in close agreement with 
preliminary results of numerical investigations, that take 
into account inductive effects 101. 
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The aim of this work is to provide insight into the ar- 
ray dynamics underlying the structures found in the I-V 
curves over the whole current range. To achieve this goal 
we compare the experimental imaging results to graphical 
animations of the time-evolution of the spatially resolved 
vortex pattern distributions obtained from our model cal- 
culations. In addition, we define and calculate a number 
of order-parameter-like quantities that characterize the 
nature of the microscopic dynamics. The main conclu- 
sion is that we can distinguish three different regions in 
the array dynamics: the subcritical region (I), the vortex 
region (II), and the region of constant differential resis- 
tance (III). In region (I) the array is in a zero- voltage 
state. Region (II) is dominated by (collective) vortex 
dynamics, contributing to the structure of the TV char- 
acteristics. In contrast, in (III) the array dynamics is 
characterized by a 'nearly coherent' behavior of the junc- 
tion phases. In Section II we describe the samples and 
the imaging technique. In Section III we introduce the 
model equations and the quantities calculated. In Section 
IV the experimental measurements and the model calcu- 
lation are discussed and compared. Our conclusions are 
presented in Section V. 



II. EXPERIMENTAL TECHNIQUES 
A. Samples 

The samples used for the present studies consist of 
two-dimensional arrays of Nb/A10 x /Nb junctions with 
square elementary cells. The junctions are square with 
an area of about 18 /jm 2 . We used 6x6, 10x10, and 
20x20 arrays without a ground plane and 10x10 and 
20x10 arrays with a superconducting Pbln ground plane 
placed at a distance to the array of about 1 /xm. Here 
an N X M array denotes an array that has N columns 
of M junctions. The lattice spacing a is 16.7 /im. Each 
of the junctions is externally shunted by an ohmic re- 
sistor R s » 1.5 Q to decrease the McCumber parame- 
ter /3 C = 27ri c i?gC/$ ~ 0.7 pj (overdamped regime), 
where $o = h/(2e) denotes the flux quantum. The criti- 
cal current of each junction is i c ~ 150 /iA. The spread 
of i c over one array is typically less than 3% (one stan- 
dard deviation from the mean value) [ fiof . The magnetic 
penetration depth = U/(2efj,oi c ) of the arrays IjllJ is 
smaller than a, where no is the permeability of free space 
and e is the elementary charge. In Ref. |l(| the sample 
geometry, layout, and fabrication is described in more 
detail. 



B. Experimental Imaging of Arrays 

Low-temperature scanning electron microscopy offers 
the possibility to image various properties of supercon- 
ducting samples during their operation at liquid helium 



The basic LTSEM princi ples together 



temperatures 

with some results are described in Refs. 



The 



top surface of the sample is scanned with the electron 
beam, while the sample is thermally coupled to a liquid 
helium bath. For the present studies, the sample is dc 
current biased and the electron beam induces a change 
AV in the array voltage that is recorded as a function of 
the focus coordinates (xo, yo) of the e-beam. In order to 
increase the sensitivity, the e-beam is chopped with a 20 
kHz frequency and AV is phase-sensitively detected with 
a lock-in amplifier. Typical values for AV are in the the 
range 100 nV - 5 fiV, whereas the array voltage, V, is of 
the order of mV. Hence, the perturbation due to the e- 
beam irradiation is small. The sample temperature is es- 
timated to be about 4.5 K for a Helium-bath temperature 
of 4.2 K. The dominant effect of the e-beam irradiation is 
local heating. We estimate from the e-beam parameters 
(electron energy: 25 keV, beam current: 100 pA) that 
there is a local temperature increment at the beam focus 
of about 0.4 K. The lateral extension of the thermally 
perturbed area is about 1 /im < a (representing the limit 
of the spatial resolution of this imaging technique). 

The sample is shielded from external magnetic fields by 
four /x-metal shields at both room and liquid helium tem- 
peratures. A perpendicular magnetic field can be applied 
using a circular copper coil placed in the liquid helium 
just below the sample substrate. In zero applied mag- 
netic field, the residual external dc field perpendicular to 
the array, corresponds to a frustration / < 0.1, with / 
the average external flux in one unit cell divided by ^o- 
This residual magnetic field B « 700 nT, was obtained 
by measuring the dependence of the array critical cur- 
rent I c on an applied perpendicular magnetic field for an 
array without a superconducting ground plane. 

The experimental data shown in this paper is obtained 
in zero applied perpendicular magnetic field, i.e. in the 
residual magnetic field mentioned above. For applied 
fields corresponding to approximately integer values of 
/ the experimental results remain qualitatively the same 
up to / w 5. 

To interpret our imaging results, it is important to 
identify the different time scales involved. The junction 
oscillation period is of the order of 10 ps, whereas the 
decay time of the beam's thermal perturbation is about 
100 ns |l2| . During step- wise scanning, the electron beam 
typically stays 3 ms at each position. The time needed to 
take one complete LTSEM image of the array dynamics is 
of the order of minutes, hence, the beam-induced voltage 
signal AV{x$, yo) represents a time- averaged quantity on 
the time scales of Josephson dynamics. 

The local temperature increment at the e-beam focus 
is most effective at the positions of the Josephson junc- 
tions. For an e-beam current of 100 pA, the heating of 
an individual junction at (xo,yo) results in a reduction 
Ai c of the critical current of this junction Ai c /i c ss 8%. 
The resulting voltage change AV(xo,yo) depends on the 
nature of the dynamics at and around the junction at 
(xo,yo)- Hence, a spatially inhomogeneous steady-state 
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dynamics gives rise to an inhomogeneous image. Spa- 
tially resolved images have been interpreted in Ref. H in 
terms of vortex motion or, in the case of underdamped 
samples, in terms of row-switched dynamical states. 

III. MODEL EQUATIONS FOR THE DYNAMICS 
OF INDUCTIVE JOSEPHSON-JUNCTION 
ARRAYS 



v = 1 V V d ^( r > r - e i/>*) ( A\ 

NMN t 2^2^ dt ■ 1 ' 

1 r t=l 

Here N t is the number of time-integration steps, ^ r is 
the sum over the N x M junctions in the current di- 
rection. V is expressed in units of i c R s , time is ex- 
pressed in units of the dimcnsionlcss characteristic time 
t c = l/u c = h/(2eR s i c ). 



In this section we briefly discuss the model we use to 
describe the array dynamics including the self-induced 
magnetic fields. In the model, the array is driven by a 
uniform applied dc current / along the vertical direction. 
In the classical regime, Ej ^> E Cl (Ej = &oi c /2ir and 
E c = e 2 /2C), the phases 6(r) of the superconducting 
order parameter on an island r, are the only variables. 
The array dynamics is then determined by the resistively 
and capacitively shunted junction (RCSJ) model for each 
junction, Kirchhoff's current conservation condition on 
the superconducting islands, plus Faraday's law for the 
magnetic field dynamics. 

Using the approximation introduced in Ref. ^] 
("Model C") for the full-range inductance matrix, and 
using the temporal gauge, one obtains a closed set of 
dynamical equations for the gauge-invariant phase dif- 
ferences *(r, r') =0(r)-0(r')-27rA(r,r'). Here A(r, r') 
is defined by the line integral of the vector potential A, 

A(r,r') = (1/<&q) J* A • dl. The derivation and imple- 
mentation of these model equations is discussed in more 
detail in Rcfs. |I|,|l4|- In the model calculations we can 
explicitly tune the dimensionless parameters (i c and Aj_, 
introduced in Section II. A. 

In our simulations we obtain ^(r, r + a) (a = e x , e y ) as 
a function of time. The magnetic flux is then given by 



2tt$(R, t) 



E *(ry,t). 



(1) 



Here 'P(R) is the anti-clockwise sum over the four bonds 
(r, r') around the plaquette with coordinate R. The vor- 
ticity n(R) is given by 



27r(n(R) - *(R)/* ) 



E *( r 

•P(R) 



r'), 



(2) 



with 



*(r,r') = tf(r,r') - 2nAf 



<f(r,r') 
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or, equivalcntly, by 
n(R, t) = 



E^ 

•P(R) 



*(r, r') 
2tt 



(3) 



Here the function Af yields the integer nearest to the ar- 
gument. The voltage response for a current applied in 
the y-direction is obtained from 



IV. EXPERIMENTAL AND THEORETICAL 
RESULTS 



A. Current- Voltage Characteristics 



Figure 1 shows a representative experimental current- 
voltage (I-V) characteristic together with its differential 
resistance dV/dl obtained from a 10x10 array. The rich 
structure of dV/dl above I c is typical for all arrays stud- 
ied. In Fig. 1 we indicate three regions. The subcritical 
region (I) defined for / < 7 C , an intermediate current 
region (II) ending at and region (III), where the 

differential resistance becomes constant, for / > Iu n . 




I(mA) 



FIG. 1. I-V characteristic of a 10x10 array together with 
the differential resistance (dV/dI)(I) measured at a temper- 
ature T » 4.5 K. 

In Fig. 2 we show an I-V curve obtained from numerical 
simulation of a 20 x 20 array with /3 C = 2iri c R 2 s C /<&q = 
and Aj_ = 0.6a. Also shown is dV/dl (thin line). The 
critical current is I c /N — 0.87i c . This I-V curve is 
qualitatively similar to the one measured in the exper- 
iments: for intermediate currents, the differential resis- 
tance shows a jagged structure, and for currents I > 
1.5Ni c the I-V curve is linear, just as is the case for 
large currents in the experimental I-V curve. As our 
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prime goal is to obtain a qualitative modeling of the ex- 
perimental systems, the choice of simulation parameters 
is motivated partly by numerical convenience. In par- 
ticular, we simulate a Ai = 0.6a, although in experi- 
ment this ratio is about 0.1. The latter value would lead 
to considerably stiffer differential equations that require 
much longer computation times. Due to their smaller 
\_i_/a ratio, the samples studied experimentally have even 
stronger self-induced magnetic fields than the simulated 
ones, and thus smaller critical currents. Furthermore, 
in experiment, region (II) extends up to Iu n ~ 2.9Ni c , 
whereas the simulations indicate that the border between 
regions (II) and (III) is at / w 1.5iVz c . 



2.0 




I/(Ni c ) 

FIG. 2. Theoretical zero temperature / — ^-characteristics 
(fat line) and the differential resistance (dV/dI)(I) (for se- 
lected data points we also show the error bars) of a 20x20 
array with f3 c = 0, and Ax = 0.6a. V is the time-averaged 
voltage across the array in the current direction. 

B. Imaging of Array Dynamics 

In this section we present images of the array dynam- 
ics for the three different regions in the I-V characteristic 
in Fig. 1. We first discuss the subcritical region, then 
the region of vortex dynamics, and finally the region of 
constant differential resistance. Typical LTSEM imaging 
results are shown in Fig. 3 and model simulation results 
in Fig. 4. 

1. Subcritical region (I < I c ) 

For a bias current I < I c and in the absence of LTSEM 
heating, the array is in the zero-voltage state. Never- 
theless we obtain useful information from the dynamical 
imaging experiments. In particular, we will see below 



that the LTSEM images obtained in this region confirm 
the importance of inductive effects in our small-A^ sam- 
ples. 

In Figure 3 (a) we show a typical LTSEM imaging re- 
sult below but close to the array critical current J c , for 
relatively high beam power. From this result we see that 
junctions at or near the edges parallel to the bias current 
give a voltage response to the local heating. This can 
be understood in the following way. For currents in this 
subcritical region, but close to the array critical current 
I c . the LTSEM is acting as an active probe, inducing 
vortex motion. When the i c of a junction at the array 
edge is lowered due to the e-beam irradiation, vortices 
can overcome the energy barrier for entry at this junc- 
tion and subsequently travel across the array. As a result, 
a voltage signal AV > is observed. The corresponding 
process for junctions not at the edges is the creation of a 
vortex-antivortex pair. From Fig. 3(a) we see that this 
latter process does not occur for junctions that are not 
close to the edges. This means that such junctions carry 
less current, which corresponds to larger energy barriers 
for vortex-antivortex creation. Such an inhomogeneous 
distribution of the bias current is in agreement to the 
small magnetic penetration depth A of the array. The 
correspondingly strong inductive effects lead to a spatial 
distribution of the dc bias current that is strongly peaked 
at the two array edges parallel to the current flow jl5],|l6| . 
The current-induced flux is also maximal at these edges. 
It is oriented in opposite direction at opposite edges. This 
is illustrated in Fig. 5, where we plot the numerically 
calculated current and magnetic field distributions for a 
20 x 20 array at I = O.S6Ni c and with Aj^ = 0.6a. Simi- 
lar simulation results in this subcritical region have been 
reported in Refs. @Q. 



2. Vortex dynamics region (I c < I < Iu n ) 

Above I c the current through the edge junctions ex- 
ceeds the junction critical current. As a result, vortices 
enter the array at one edge and antivortices at the oppo- 
site edge. These vortices are depinned from the edges by 
the Lorentz force and move across the array, generating 
the observed voltage. 

LTSEM results for bias currents slightly above I c are 
presented and discussed in Refs. In the LTSEM 

images obtained in this current region it is seen (see, 
e.g., Ref. ||) that the sign of the voltage signal near the 
sample edges tends to alternate along the current direc- 
tion. These images indicate an alternating or staggered 
crossing vortex motion where vortices and antivortices 
are nucleated at opposite array edges and subsequently 
move across the whole array. At the array edge opposite 
to the nucleation site they leave the array or, equivalently, 
annihilate with an image vortex of opposite sign. 
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(a) (b) 




x (|um) 150 Xo(nm) 

FIG. 3. Grey value representations of the experimental voltage image AV(xo,yo) for a 10x10 array at T ~ 4.5 K. The 
array is current biased at I = 0.9J C (a) and / = 4.25/ c (b), respectively. The dc bias current flows vertically through the 
array. The array boundaries lie between /jm and 150 fim in both directions. A positive (negative) e-beam induced voltage 
signal AV(xo,yo) is indicated by the dark (bright) areas, whereas zero signal is shown by the area surrounding the array. The 
individual rows of junctions are indicated by the small arrows numbered 1-10 from top to bottom. In (a) , the A V < voltage 
response at the top and bottom of the array columns (marked by arrows from the left) arises from the current feeding resistors 
made from InAu thin films. 



(a) (b) (c) 




FIG. 4. Vortex configurations for simulations of a 20x20 array, for three different values of the bias current (a-c). For each 
current value, we show four consecutive frames (0-3). The dc bias current flows vertically through the array. The black (white) 
squares denote plaquettes with vortex number n(R) = +1(— 1). Vortices with positive sign move to the right, vortices with 
negative sign to the left, as can be deduced by comparing consecutive frames, (a) I = 0.95iVi c . Frame is at t = 3000t c . 
Between consecutive frames there is a time interval of 3i c . Adjacent rows tend to be crossed by vortices of opposite sign. The 
alternating structure is disrupted between the third and fourth row from above and between the sixth and seventh row from 
below, (b) I — 1.20iVi c . Frame is at t = 9025t c . Between consecutive frames there is a time interval of 2t c . (c) I — 2.0Ni c . 
Frame is at t = 500.75t c . Between consecutive frames there is a time interval of 0.9i c . 
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We will see below that this interpretation is confirmed 
by the general behavior found in simulation images for 
a bias current slightly above the array critical current. 
This kind of alternating crossing vortex motion has no 
observed analog in continuous superconducting samples. 

To analyze the dynamics of the vortices in full detail, 
we have studied graphical animations of the time evolu- 
tion of the vortex distributions in the array. First we dis- 
cuss the ones obtained in a simulation for / = 0.95iVi c of 
a 20 x 20 array with j3 c = and = 0.6a. When we start 
the simulation with random initial phases we observe, af- 
ter a transient of about t/t c — 800 time units, the type 
of vortex motion depicted in Fig. 4 (a). The snapshots 
in Fig. 4 (a) show the same type of vortex motion as the 
one deduced from the LTSEM measurements: vortices 
of opposite sign cross the array in opposite directions 
in adjacent rows. We observe that the staggered struc- 
ture is broken at two places, where two adjacent rows 
are crossed by vortices of equal sign. One might view 
these places as domain wall defects between two differ- 
ent polarities of the staggered pattern (for experimental 
results see also Ref. ||). We have verified that this type 
of dynamics is stable for very long simulation times. The 
position and number of domain walls depends sensitively 
on the initial conditions. For all currents in the range 
I c /(Ni c ) = 0.87 < I/(Ni c ) < 1.15 the long-time stable 
vortex dynamics is of the same staggered type. 




R 

FIG. 5. Simulated field and current distribution in the cen- 
tral row of a 20 x 20 array with. Ax = 0.6a and / = 0.86iVi c , 
just below the array critical current. i y is the current through 
a longitudinal junction. 

The dynamical patterns observed experimentally as 
well as numerically for / > I c are strongly influenced 
by vortex- vortex interaction. These interactions lead to 
the almost regular patterns in which the vortices tend to 
move. Towards the high current end of region (II), the 
dynamics is different. Imaging the dynamical state by 
LTSEM now yields larger two-dimensional domain pat- 
terns spreading over several unit cells in both the x and 
y-directions. A typical example for this kind of voltage 



response is given in Fig. 3 (b). From our experimental 
observations, we can deduce the following: (1) The dy- 
namical state, which gives rise to the voltage image is 
stable in time. If the parameters of the sample (I, T) 
are not changed a subsequent image will give the same 
result. (2) The imaging results depend sensitivily on the 
history. If, for example, the bias current or tempera- 
ture is changed significantly, and then returned to the 
same parameter values, the imaging results change. (3) 
When / or T is changed smoothly, we observe a smooth 
variation of the detected patterns. (4) A magnetic field 
(noninteger / or |/| > 5) changes the regular pattern, ob- 
served for small bias currents discussed above, to a com- 
plex response similar to that shown in Fig. 3 (b). Based 
on these observations, we conclude that the voltage re- 
sponse is caused by a complex multi- vortex dynamics, 
and not by e.g. sample inhomogeneities, trapped flux, or 
temperature fluctuations. 

In our simulations, the staggered vortex dynamics is 
the relevant dynamics up to approximately / = 1.15iVi G . 
The region 1.15 < I/(Ni c ) < 1.5 is a transition region 
between the regime of staggered vortex dynamics and the 
regime of constant differential resistance. In this current 
interval, in some parts of the array vortices move inde- 
pendently, and in others we observe vortices that tend 
to move coherently in adjacent rows. A representative 
example of this type of dynamics is given in Fig. 4 (b). 

3. Linear branch (I > Iu„) 

In region (III), where the I-V curve is linear (con- 
stant differential resistance), each junction in the array 
columns (longitudinal junction) yields approximately the 
same voltage signal AV(io, yo) and the LTSEM image is 
rather uniform over the whole array. For example, at 
/ = 5I C for a 10 x 10 array the voltage signal at each 
junction was the same within 5% (see Fig. 7 of ||). 

The absence of any structure in the beam-induced volt- 
age signal suggests that there are no isolated vortices en- 
tering or leaving the array. This is indeed what we find 
in the simulations in the region of constant differential 
resistance. In fact, the dynamics in this current region is 
not due to vortex motion but to a state in which the lon- 
gitudinal junctions belonging to the same column oscil- 
late almost in phase. This phase coherence reveals itself 
in a wave-like dynamics of the magnetic field distribu- 
tion. When looking at the discrete vortex configurations 
shown in Fig. 4 (c), we observe fronts of vortices that 
move inward from the boundaries. In the middle the 
vortices annihilate. The discrete vortex configurations in 
Fig. 4 (c) show a high degree of symmetry, which reflects 
the coherence in the motion of the longitudinal junctions 
in different rows. We have also explicitly simulated a 
scanning-induced 8% critical-current reduction of subse- 
quent individual junctions for / = 2.0Ni c , and indeed 
find a spatially uniform voltage change. 
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FIG. 6. Crossover from vortex-dominated to coherent-phase dynamics as a function of applied current /, probed by (a) V 
(full line), V/n v (dot-dashed) and the variance an of the Hall voltage (inset, in arbitrary units); (b) P (full line) and P s 
(dot-dashed) and n v (inset). All results are obtained in an upward current sweep of a 20 x 20 array for k = 0.6 and T = 0. For 
each current value we used a warm-up time of 1000t c and an averaging time of 2000t c . 



C. Crossover from vortex-dominated to 
coherent-phase dynamics 

In the above we have focused on comparing the LTSEM 
images with snapshots of vortex configurations from nu- 
merical simulations. We found three regions in the I —V 
characteristics as well as in the experimental and numer- 
ically obtained images. We related these regions to dif- 
ferent types of dynamics. The crossover from a vortex- 
dominated to a coherent-phase dynamics can be explored 
in more detail in the model simulations. To this end, we 
can consider a number of quantities that probe the degree 
of vortex organisation or the degree of phase coherence. 
In particular, we find that the following order-parameter- 
like quantities P and P s can be used to distinguish the 
different current ranges, 



Ey=i Mr,*) I 

N v {t) 



Ey=i(-i) y n e (r,*) ' 

N v (t) 



(5) 



(6) 



where 



JV-l 



n c (Y,t)= ^n(R,t), N v (i) = ^ |n(R, t) \ . 



R 



X and Y denote the x and y component of the plaqucttc 
coordinate R, respectively: R = Xe x + Ye y . The phys- 
ical meaning of P and P s can be inferred from the fact 



that a vortex that enters at one side of the array, either 
leaves the array on the other side, or is annihilated by an 
antivortex moving in the opposite direction. The quan- 
tity n c (Y,t) for row Y distinguishes between these two 
possibilities. In the former case n c is nonzero, whereas in 
the latter case it is, on average, zero. A value of P = 1 
thus implies that all the vortices cross the whole array 
unobstructed, while P = implies that the vortices are 
annihilated in the middle of the array. The staggered or- 
der parameter P s measures whether the spatially resolved 
vortex dynamics consists of alternating rows of vortices 
and antivorticcs crossing the array. The presence of do- 
main walls, i.e. two adjacent rows in which vortices cross 
the array in the same direction, reduces P s . 

On the basis of this interpretation we expect that Fig. 
4(a) (low current region (II)) corresponds to P w 1, 
P s w 1; and Fig. 4(c) (region (III)) to P « 0, and P s w 0. 
In Fig. 6 we plot P, and P s versus the applied current. In- 
deed, directly after depinning, J > I c « 0.87Ni c , P and 
P s attain values close to one. As the current increases 
this value slowly decreases (low current region (II)). At 
I w \.2Ni c the value of both order parameters exhibits a 
sharp drop to a much smaller (but nonzero) value. In this 
region the order parameters slowly decrease to the values 
near zero in (III). Thus, using P and P s one can readily 
establish the type of vortex dynamics without having to 
study vortex animations for each value of the bias cur- 
rent. In particular, long-time-stable values P w 1, P s w 1 
correspond to a staggered crossing vortex motion. 

We will now correlate the different regions, mapped out 
using P and P s , with the behavior of other quantities. In 
Fig. 6(a) we show both the voltage and the average ve- 



7 



locity per vortex V v = V/n v (the voltage normalized by 
the number of vortices) versus current. In the inset of 
Fig. 6(b) we plot the vortex density. We observe that 
the vortex density displays a pronounced maximum in 
region (II) , accompanied by a smaller slope of the V v (I) 
curve as compared to the other regions. The dynamical 
properties of the vortices in this region are thus different 
from the other regions. 

Vortex jumps across a junction do not only give rise to 
a contribution to the longitudinal voltage across the ar- 
ray, but also induce a fluctuation of the transverse (Hall) 
voltage around zero. For the nearly antisymmetric dy- 
namical patterns of vortices in the phase-coherence re- 
gion (III), like the ones shown in Fig. 4(c), the Hall volt- 
age contributions of the left and right half of the array 
tend to cancel by symmetry. In the region with stag- 
gered crossing vortex motion however, the vortex jumps 
in the left half of the array do not occur in unison with 
antivortex jumps in the right half, leading to larger fluc- 
tuations in the Hall voltage. Therefore the Hall voltage 
fluctuations may be viewed as a measure of the degree of 
(anti)symmctry in the dynamics. Indeed, as seen in the 
inset of Fig. 6(a), the magnitude of these fluctuations is 
reduced dramatically between the vortex-dominated and 
the coherent-phase regime. 

Figure 6(a) is obtained in an upward current sweep. 
We have also performed a downward current sweep from 
region (III), starting with a uniform phase configura- 
tion. We find that region (III) is exactly characterized 
by P = P s =0, whereas in the upward sweep there are 
still some asymmetries in the vortex configuration that 
yield small but nonzero values for the order parameters. 
Entering the high current end of (II), P and P s attain 
nonzero values. There is a slight hysteresis in the cur- 
rent value at which P and P s attain values close to one. 
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The power spectrum of the voltage, 
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FIG. 7. Calculated power spectrum of the array voltage for 
three different bias currents. Subsequent curves are offset by 
5.5 units. 



dt V(t)e 



%2-KUt 



(7) 



is another useful probe for the vortex dynamics. A non- 
zero voltage can be viewed as either created by a vor- 
tex jump, or by a junction phase slip. When the lon- 
gitudinal junctions oscillate coherently, S{y) consists of 
sharp peaks at multiples of v = V/ (2irt c ). If the dynam- 
ics, however, is dominated by incoherent vortex jumps, 
the peaks in the spectrum are much broader. In Fig. 7 
we show S(y) for 3 different currents. From if, = 0.95, 
within the "vortex" region, to i = 2.0 (region (III)) the 
spectrum changes from a noisy to sharply peaked one. 
These results again illustrate the crossover from a vortex- 
dominated dynamics to a coherent-phase dynamics for 
higher currents. 



V. CONCLUSIONS AND DISCUSSION 

In summary, we have shown that the comparison of LT- 
SEM images and results of model calculations can signifi- 
cantly increase our insight in the dynamics of Josephson- 
junction arrays. On the one hand, the comparison of the 
experimental images with the results of model simula- 
tions corroborates the interpretation of images for cur- 
rents not too far from the array critical current, and con- 
tributes to our understanding of the dynamics underly- 
ing the images outside this current region. It shows that 
the dynamics we found exists even in the absence of the 
disturbance produced by the measuring device. On the 
other hand, the agreement with the experimental results 
supports the relevance of the model equations employed. 
In the context of the samples studied here, an essential 
ingredient of these model equations is the inclusion of 
(strong) self- induced magnetic fields. 

The successful comparison between experimental and 
theoretical results has enabled us to map out three re- 
gions with different types of dynamics for dc biased arrays 
with Aj^ < a in zero applied magnetic field. For bias cur- 
rents above the array critical current we found current- 
induced vortex nucleation at the array edges parallel to 
the bias current. We identified an alternating pattern of 
crossing vortices and antivortices as the typical vortex 
dynamics existing for bias currents slightly above the ar- 
ray critical current. We conclude that at least part of 
the rich structure found in the experimental J-V's is due 
to the dynamics of vortices ]i7|JT8| . For larger currents, 
the I — V characteristic becomes linear, and the under- 
lying dynamics is characterized by a growing tendency 
of longitudinal junctions to oscillate in phase. We have 
further illustrated the crossover from vortex-dominated 
to coherent-phase like dynamics by numerically studying 
the spectral function, the Hall-voltage fluctuations and 
order-parameter-like quantities. 



8 



Recently Oppenlander and coworkers jl9[ have also cal- 
culated I-V characteristics similar to the experimental 
ones, however without taking into account the dynamics 
of the self-induced magnetic fields that are essential for 
a successful theoretical description of the experimentally 
observed microscopic dynamics. In their work the struc- 
ture in the differential resistance is due to the inclusion 
of strong (magnetic) disorder [2(J ■ 

In this paper, we have discussed the numerical results 
for zero McCumber parameter. In the experiments, the 
McCumber parameter was estimated to be (3 C = 0.7. In 
our model simulations, the microscopic dynamics for cur- 
rents above but close to the array critical current remains 
qualitatively the same for this value of the McCumber 
parameter. For higher currents the dynamics is again 
characterized by a large degree of spatial coherence. The 
region of genuine vortex dynamics shrinks with increasing 
(3 C . For (3 C > 2.5 the system enters a row-switched state 
immediately above the critical current (for Ax = 0.6a), 
and therefore no vortex-flow regime is found. For non- 
zero applied magnetic field, row-switched states in induc- 
tive arrays were studied in Ref . . 

For the arrays used in the present studies, the 
Kosterlitz-Thouless-Berezinskii phase transition temper- 
ature Tktb is close to the superconducting transition 
temperature T c of the Nb thin films. The experiments are 
performed at temperatures well below Tktb, 2e ^ B T ~ 
1 • 10~ 3 . For these temperatures the effect of thermally 
induced vortices on the array dynamics is negligible p2[ . 
We have also performed calculations for \± > a (outside 
the region of the present experiments) up to X± = 10a 
p3| . In this regime we have found a staggered vortex 
dynamics similar to the one observed for X± ^ a, again 
occuring for currents slightly above the array critical cur- 
rent. This indicates that the staggered vortex dynamics 
is the generic dynamics for such currents. 

In the LTSEM experiments, the images correspond- 
ing to the staggered vortex dynamics were found for 
some range of magnetic frustrations around zero. We 
have simulated the dynamics in the lower part of re- 
gion (II) (where the alternating crossing vortex motion 
is observed) in a magnetic frustration of / = 0.01. For 
currents slightly above I c , we find that the array evolves 
towards a similar state as found for / = 0, but now there 
are on average more vortices with positive vorticity than 
with negative vorticity. 
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